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ADAPTIVE CONTROL OF LARGE SPACE STRUCTURES 
USING RECURSIVE LATTICE FILTERS 


By 

Gene L. Goglia* 

SUMMARY 

This report summarizes the research activities performed under grant 
NAG-1-429. The objective of the research has been to study the use of re- 
cursive lattice filters for identification and adaptive control of large 
space structures. Lattice filters are used widely in the areas of speech 
and signal processing. Herein, they are used to identify the structural s 

dynamics model of the flexible structures. This identified model is then 
used for adaptive control. Before the identified model and control laws are 
integrated, the identified model is passed through a series of validation 

procedures and only when the model passes these validation procedures con- >■; 

U 

trol is engaged. This type of validation scheme prevents instability when ij 

the overall loop is closed. !j 

One of the main aims of the research has been to compare the results 

[! 

obtained from simulation to those obtained from experiments. In this re- 
gard, the flexible beam and grid apparatus at the Aerospace Control Research 
Lab (ACRL) of NASA Langley Research Center were used as the principal candi- 
dates for carrying out the above tasks. Another important area of research 
namely that of robust controller synthesis was investigated using frequency 
domain multivariable controller synthesis methods. The method uses the 
Linear Quadratic Gaussi an/Lobp Transfer Recovery (LQG/LTR) approach to en- 
sure stability against unmodeled higher frequency modes and achieves the jj 

*Eminent Professor, Department of Mechanical Engineering and Mechanics, Old 
Dominion University, Norfolk, Virginia 23508. h 


desired performance. Such a controller was designed for the 122 mr. Hoop- 
Column antenna using a single 3-axis torque actuator and attitude sensors. 

References 1 and 2 present the detailed analysis of identification 
results for the flexible grid apparatus using lattice filters. The scheme 
provides on-line identification of nunber of modes, mode shapes, modal damp- 
ing and natural frequencies. The results indicate that the lattice identi- 
fication scheme is a viable scheme for identifying the structural dynamics 
of flexible structures. The experimental results also indicate differences 
between those predicted by finite element analysis and obtained by experi- 

t 

ments. The difficulties are not as such in finite element analysis but in 
modeling the apparatus for finite element analysis. This fact emphasises 
the need for on-orbit identification of large space structures before con- 
trol is attempted. A summary of the experimental results obtained using 
lattice filters is described in reference 3. 

An adaptive control scheme using lattice filter identification and 
modal description has bean developed in reference 4, Alternate schemes of 
using input-output models instead of modal form from lattice filters is 
described therein. The problem in this approach is to obtain efficient 
control schemes as the identified model of the system becomes coupled arid to 
calculate the pole placement control law on-line is computationally complex. 
Presently, the identification scheme using lattice filters for obtaining the 
input-output model is under development in the Charles River Data Systems in 
the ACRL. 

A new approach of designing robust controller for a large flexible 
space antenna using the LQG/LTR approach was developed in reference 5. The 
method was used in designing robust controller for the 122 mr. Hoop-Column 
antenna using only a 3 axis torque actuator and attitude sensor. The objec- 


2 


tive is to design the controller based on a lower order model to achieve the 
desired bandwidth and at the same time ensuring stability against unmodeled 
higher frequency modes. The results In reference 5 indicated that If one 
uses only a rigid body model for design stability against unmodeled modes 
can be obtained but not the performance. Based on detailed studies in ref- 
erences 6 and 7 if was concluded that with the first three flexible modes 
(corresponding to the 3 axes) included in the design model both stability 
and performance can be ensured. 
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Abstract 

The basic theory of least squares lattice filters and their use in identification of structural 
dynamics systems is summarized. Thereafter, this theory is applied to a two-dimensional grid 
structure made of overlapping bars. Previously, this theory has been applied to an integral beam 
System identification results are presented for both simulated and experimental tests und they arc 
compared with those predicted by means of finite clem ..it modeling, The lattice filtering ap- 
proach works well for simulated data based on finite element modeling. However, considerable 
discrepancy exists between estimates obtained from experimental data and the finite element 
analysis, It is believed that this discrepancy is the result of inadequacies in the finite clement 
modeling to represent the damped motion of the laboratory' apparatus, 

Introduction 

The ability to predict the dynamic behavior of large space structures (LSS) ade- 
quately for control system design is doubtful because of their expected size, appreciable 
flexibility, and on-orbit assembly anomalies, Hence, dynamical modeling from on- 
orbit measurements, followed by modifying the control system as dictated by the 
identified control system design model (adaptive control), is of interest, The goal of this 
paper is to determine, using a generic grid structure, whether a priori modeling of the 
structure is adequate for a high authority control system design or whether on-orbit 
identification is needed. 

An approach for identifying the dynamic behavior of LSS that estimates model order 
in addition to model parameters is presented in [I). It uses lattice filters which provide 
an order as well as a time recursive algorithm for linear least squares signal estimation. 
{2] provides a tutorial on lattice filter theory and applications. The outputs of the theory 
of [1] are the least ‘square estimate of the measurement sequence, the model order 
required to fit the measurements, the associated lattice model (this includes mode shape 

'This paper was originally presented ai the American Control Conference. San Diego, California, June 6-8. 
1984. 

'NASA Langley Research Center. Hampton, VA 23665, 

'Old Dominion University Research Foundation, Hampton, VA 23666 
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estimates that are onhonormal in the measurement space), and the associated auto* 
regressive moving average lARMAX) model of the measurement sequence. 

Some distributed adaptive control strategies require identification of the natural 
modes of a structure |3. 4|. Unfortunately, the lattice filter provides mode shapes that 
are orthonormal in the measurement space and, hence, are not the natural modes. 
Natural modes can be obtained, however, either through an eigenvector analysis of the 
identified ARMAX model or through a transformation that provides spectral decom* 
position of the lattice filter modal amplitudes (5). The latter method is nsed herein to 
obtain the natural modes. Using spectral decoupling to determine mode shapes, one can 
obtain mode frequency and damping using an equation error parameter identification 
method [3] that employs a second-order ARMAX model to represent the natural mode 
amplitudes, The procedure of (3| tracks frequency and damping coefficients required 
for the modul amplitude input sequence to fit the second order ARMAX model. 
Because the input sequence may have an unfavorable signal to noise ratio, the parame- 
ters so derived must be tested for validity before use in control system design. This 
parameter testing is treated in detail in (6| but is not employed herein. 

The foregoing procedure has been used to identify the dynamic characteristics of an 
integral free-free beam in [l.7|. |8| describes the test apparatus used in those studies. 
In this paper, the theory is applied to a more complex, two-dimensional grid structure 
made of overlapping bars, First, a brief overview of the theory used is presented. Next, 
system identification results are presented using both simulated and experimental 
data. Finally, the experimental results are compared, with those predicted using finite 
element modeling, 

Summary of the Method Used to Identify Structural Dynamics Systems 

For the application considered here we assume that the A’th measurement sample is 
of the form 

yl ~ y\yU')J (I) 

where NS represents the number of sensors. It is assumed that y is generated from a 
model system such that 

y* “ + tn (7.) 

Here, d> is a mode shape matrix, 'Ft is the modal amplitude vector, und m is a Gaussian 
random variable with a zero mean and a covariance matrix R, Typieully. for structural 
dynamics applications, each component of 'R* is the output of an uncoupled second- 
order process. The task here is to estimate the order and obtain the least square estimate 
of from N + I measurement samples v 0 through v\ , [ 1 ] presents a derivation of the 
equations that relate order n, and time i. recursions for the normalized forward and 
backward residuals ( e and r, respectively) as well as the least squares estimate of the 
measurement vector y. These equations are listed below: 

<Y«*i = (1 ~ Arr.n-i)* 1 J (e, „ - (3) 

| f l f — | . /, ” 

,V-I 

y.v = £ E(e.v.,!r lV -u) 

n“0 


( 5 ) 
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with 

*1 •. I " (e M .r, -!.»,> (6) 

and B{x i.v) is the orthogonal projection operator of the vector .v onto the vector y. The 
symmetry of the recursion formulae is apparent. The equations are coupled by the term 
A’, „,| which is customarily called the “reflection coefficient." 

Clearly, in tins approach one may ’fit the noise" by continually increasing the order 
of the system; however, once the order of the estimator has increased sufficiently, the 
residual errors should lie within a noise band which can be predicted based on assumed 
noise characteristics. A threshold value can be selected based on this predicted noise 
band and order determined by a test of whether or not the residuals have been reduced 
to lie within the noise band. Also, the test can be made considering several samples of 
data; that is. using a data window. [I] documents experience in order determination 
based on this threshold test. 

Having defined the order required to fit the data using a linear model, we seek a fixed 
set of basis functions that are spectrally decoupled for modal control, Therefore, u fixed 
orthonormal basis is used during intervals when the order estimate is constant. (How* 
ever, the order estimate is checked at each measurement sample based on the threshold 
test.) The lattice filter uses the current measurement as the first mode shape and, using 
a modified Gram*Schmidt orthonormalization procedure, generates additional basis 
functions from estimation residuals. Consequently, the output of the lattice filter pro* 
duces coupled mode shapes and corresponding modal amplitudes wherein the first 
coupled modal amplitude will contain ojl significant natural modes. Since the order 
estimate n has been determined, the first coupled mode digital Fourier transform (DFT) 
amplitude spectrum is searched for the n most significant peaks and corresponding 
frequencies. Because the spectrum contains n peaks for the n separate modes, a 
transformation matrix can be obtained that decouples the spectrum. This transformation 
matrix is the inverse of the matrix whose elements are the real pan of the transform of 
the n coupled modal amplitude channels (rows) evaluated at the n peak frequencies 
(columns). It effectively transforms the lattice filter modes into spectrally decoupled 
natural modes. These decoupled modes are not onhogonal. This procedure is described 
in [5], 

Thus, the decoupled modal amplitude time scries, w(k), is obtained by applying the 
transformation to the direct output of the lattice filter. I his time series is then analyzed, 
for each mede, to identify the parameters of its autoregressive moving average 
ARMAX model. The inputs to each ARMAX modal model arc the generalized forces 
and hence, each model takes on the form 

w(k) * A,w(A- - 1> + A 2 h# - 2) + BJ(k - I) + B : f(k - 2) (7) 

where/ represents generalized forces. The parametersp r = (Ai.Aj.fli.#;) are the ones 
which are identified and which are required for the control law design process. Thus, 
the ARMAX model output error is 

e(k - 1) - wlk - 1) - [A i (A* - IMA- - 2) + A : (A- - IMA- - 2) 

+ £,<A - 1)/(A- - 2) + Bi(k - \)f(k - 3)] (8) 

The method of [3] is used to identify the parameters (p) using the iteration sequence 
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p(k) * p(k - 1) + e{k - 1) 

* [WMk - 2),WMk - 3), W)/(k ~ 2 ),W t f(k - 3)) (9) 

As indicated in [3], the weights W, (l ■ 1 ..... 4) must be selected so that they arc 
consistent with the relation 

- 2) + Wiw'ik - 3) + W } fHk - 2) + W t ,/ 2 (k - 3) < 2 (JO) 

and the inputs to the algorithm (w and/) must be sufficiently varying and large if the 
parameters arc to converge to their correct value. 

Damping ratios and natural frequencies of the modes can be obtained from A i and 
Ai of equation (7), However, the solution is not unique, due to the foldovcr phenome* 
non of sampling. By finding the roots, z„ of the characteristic equation (7) and using 
the relation z, “ e' 1 ’ in the primary strip, where r is the sampling period, the following 
equations for damping ratio and natural frequency arc obtained for a typical root,* 
say z; 


u) * <j)/(2nr) 

(ID 

C « c(c* + d» J ) l,J 

(11. 

d> « tan ’Kb/a) 

(13) 

c * ln(o J + b : ) 

(14) 


and a * Re(z) and b ■ Im(z). The behavior of this overall system identifica* 
tion methodology with both simulated and experimental data is discussed in the sub- 
sequent sections. 

Description of the Flexible Grid Facility 

Figure l shows the flexible grid experimental appartus currently being built at the 
NASA Langley Research Center. The grid is a 7 ft by 10 ft planar structure made by 
overlaying aluminum bars of rectangular cross section. The bars are centered every foot 
so that there are 8 vertical and 11 horizontal bars, As shown in Fig. E, the grid is 
suspended by a cable at two locations on the top horizontal bar. The motions of the grid 
perpendicular to the plane of Fig. 2 are the ones of interest in this study. There are nine 
noncontacting deflection sensors mounted on a back frame which give a 9 x 1 
measurement vector. The sensor data are linked to the main CYBER 175 Real-Time 
Computer System at NASA Langley Research Center so that the identification can be 
carried out in real time. For the experimental tests, the locations of the sensors are 
indicated in Fig. 2. 

Simulation Studies 

A finite element analysis of the grid was performed which included the suspension 
cables. Nodes W'ere placed at each overlapping joint on the grid, the ceiling attachment 
points of the cable, and every one-half foot along the cable. The grid elements con- 
necting the nodes were modeled as bending elements, whereas the cable elements were 
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FIG. 2. A Schematic of the Grid Apparatus Indicating Locations Referred to in the Text and 

Subsequent Figures 


modeled as two-force members. Thus, a total of 165 elements were included in the 
model, Four degrees of freedom appropriate for motion normal to the plane of the grid 
were considered. No damping was included in the model. Thirty modes were obtained 
from this analysis. The frequencies of the first ten modes are listed in Table I . The first 
three modes are the pendulum modes, the fourth is the first bending mode and the fifth 
is the first torsional mode. The finite element analysis uses an iterative method to 
calculate mode frequencies. The frequencies used in simulation are believed to be 
numerically accurate since the change in eigenvalue iterate of the highest frequency 
mode used in the simulation is 10' 11 on the final iteration. The corresponding eigen- 
value iterate was 1565, which corresponds to the mode 8 frequency in Table I. 

A simulation was developed that accommodates the first 15 modes of the analysis, 
but only four m©u\?s were used herein. Modes 4, 6. 7, and 8 were used. A sampling 
rate of 32 Hz was simulated with a standard deviation for the measurement noise of 
0.005 in. which was. based on actual sensor characteristics, Modes were simulated with 
modal amplitude initial conditions of 0.1. The data window for order determination 
included eight samples. In this work, the sensor locations were chosen based on several 
simulations. These locations differ from those of the experimental apparatus in that they 
were selected to maximize the ' Tfect of simulated modes on the sensors. This was 
accomplished by visual examination of the simulated sensor outputs. The selected 
locations are indicated in Fig, 2. An asterisk is used to distinguish simulation sensor 
locations from experimental ones. One may expect that location 5 would be preferable 
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TABLE I. Modal Frequencies Obtained 
from the Finite Element Analysis of 
the Grid 


Mode Number 

Frequency (Hz) 

1 

0 364 

2 

0.625 

3 

1.398 

4 

2,29 

jt 

3.07 

6 

4.791 

7 

5 933 

8 

6.297 

9 

7.337 

10 

10.352 


to location 5**, but, since some simulated modes had little input to a sensor at loca- 
tion 5, location 5* proved to be ,5 better location. 

Based on the entire measurement vector, the lattice filler order estimate is shown in 
Fig. 3. Also. sensorS* data, typical of those of the other sensors, is shown in Fig. 3. 
After estimating the order, we carried out a transformation based on the discrete Fourier 



0 2 4 6 8 10 

TIME, SEC. 

FIG. 3. Simulation Time Histories of Sensor 5*. Lattice Filter Order Estimate, and the Norm of the 
Estimation Error for the Entire Measurement Vector. 
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FIG. 4. Sensors I and 4 from Experimental Tape 5 and the Estimate of Signal Order Obtained by the 

Lattice Filter. 

transform (DFT) using 128 samples in order to obtain the natural modes, and used the 
equation error method to identify associated modal frequencies and damping. The 
resulting modal frequencies, damping and mode shapes are compared with those 
predicted by finite element analysis in Table 2. The identification of frequencies and 
damping are close for all four simulated modes. However, the mode shape estimates 
agree with simulation for only three modes. One possible explanation for this is the 
limitation imposed by sampling rate and the number of samples used to decouple the 
lattice filter modes. Sampling at 32 Hz and including 128 data points in the DFT, a 

TABLE 2. Comparison Between Simulated and Identified Results 
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frequency resolution of only 0.25 Hz is obtained. Since the expected frequency sepa- 
ration between modes 7 and 8 is only 0.4 Hz, good decoupling cannot be achieved. 

To summarize the lessons learned from the simulation studies: 

1. The least squares lattice filter gives good identification of simulated modal 
frequencies, damping ratios, and mode shapes in the presence of sensor noise 
expected in the experimental apparatus. 

2. The DFT method of obtaining natural modes from the lattice modes is inaccurate 
if the modes are closely spaced in frequency. This may be improved by adding 
more samples to the DFT, 

3. Sensor locations should be properly selected to insure good identification of 
simulated mode shapes. 

The next section will discuss results obtained from the experimental apparatus. 

Experimental Results 

Experiments were conducted using the grid apparatus previously described, The grid 
was excited using an air shaker which periodically exhausted a jet of air that impinged 
on the grid at sensor location 1 , The frequency of the jet was adjustable from 0 to 
50 Hz. Although the resulting grid excitation was periodic, it was not purely sinusoidal 
hut was rich in harmonics, Because of the range limits of the deflection sensors - 0 to 
approximately 2 in, - the maximum peak-to-peak deflections of the grid were limited 
to about 1 in, When the peak-to-peak deflection neared this limit, the air shaker was 
turned off and the grid was allowed to vibrate freely with only air and material damping. 
A CYBER 175 Real-Time Computer System sampled the deflection sensor data at 
32 Hz for 5 seconds. The data were stored on a system data file for further analysis. 
Since only free-decay response data were recorded, the and B : parameters of 
equation (7) were not identifiable, Figure 4 presents data from file 5. Here, the order 
estimate is seen to converge to an oscillation between 2 and 3 at about 0,5 s. At about 
0,8 s, the order estimate was fixed at 3 and data collection (at 32 Hz) for the 64 time 
samples required for the DFT was started. The DFT was accomplished at about 2.8 s 
and the decoupling transformation matrix was calculated. The modal amplitudes after 
this time should contain a single frequency and the transformed mode shapes should 
correspond to the excited natural modes of the structure. Thus, three modes were 
extracted from the experimental data tape. These have frequencies near 0.5 Hz. 
2.5 Hz, and 5 Hz. Table 3 presents the mode shape estimates obtained from the 
experiment, Also presented are selected mode shape predictions taken from finite 
element analyses. The modes selected were those whose frequencies bracket the experi- 
mentally derived ones, The following discussion deals with the Table 3 data in order 
of increasing frequency, 

A good comparison does not exist between the first or third experimental modes and 
either of their bracketing finite element analysis modes. Additionally, there is some 
bending in the first experimental mode as evidenced by sensors 4, 5, and 6. The 
amplitude of this mode is shown in figure 5 along with its ARM AX parameters A, and 
A : and their primary strip equivalents of damping and frequency. Figures 6 and 7 show 
the same information for the second and third modes, respectively. For the second 
mode, a good comparison does exist between it and the 3.07 Hz finite element analysis 
mode. Note that, however, the output of sensor 4 is opposite in sign and reduced in 
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FIG, 6. Characteristics of the Second Mode Identified from Experimental Tape 5. 


amplitude from the finite clement prediction. This means that a feedback on that sensor 
based on the finite element analysis will be destabilizing near the 2,5 Hz frequency, 
The validity of this deduction can be established by examining the outputs of sensors 1 
and 4 (4). According to the finite clement analysis, the 2.5 Hz content of the sensors 
should be opposite in sign. However, they are in phase in agreement with the identi- 
fication results. 

Conclusion 

The application of the least squares lattice filter in system identification has been 
extended to a non-integral, two-dimensional grid structure made of overlapping bars. 
Previous experience has been limited to an integral free-free beam. Both simulation and 
experimental data were used to evaluate the system identification capabilities of the 
method, In the simulations, the least squares lattice filter gave good identification of 
simulated modal frequencies, damping, and mode shapes in the presence of sensor 
noise expected in the experimental apparatus. However, the spectral decoupling method 
of obtaining natural modes from lattice filter modes required a large number of data 
points in the discrete Fourier transform to get adequate frequency resolution when the 
modal frequencies were closely spaced. This problem can be overcome by an eigen- 
vector analysis of the lattice filter’s associated ARM AX model. When the lattice filter 
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FIG. 7. Characteristics of the Third Mode Identified from Experimental Tape 5 


was used for system identification with experimental data, the mode shapes identified 
differed significantly from those of the finite element analysis. This has been corrobo- 
rated by examination of the sensor data and indicates that on-line identification of large 
structural dynamic systems may be absolutely necessary to get acceptable performance 
in a high gain system that requires knowledge of mode shapes. 
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Abstract, This paper presents the use of least square lattice filters In 
identification of the dynamics of highly flexible structures. Lattice filters have 

been used extensively In the areas of adaptive signal proceslng and speech synthesis. 
Herein, they are used for on-line Identification of the nuaber of oodes, node shapes 
and aodal aaplltude tine series fro* the aeaaureaent data. The theory la lllustratad 
using experiaental data for a staple free-free beaa and a sore coaplex, flexible, 
tvo-diaenslonal grid apparatus. Results presented Indicate that tlte lattice filter 
approach produces effective Identification of structural dynamics for the class of 
structures studied to this tine. 

Keywords t Identification, Lattice Filters, Space Vehicles, Least Squares 

Approxlaation, Vibration Measurements, 


INTRODUCTION 

With the site of the structures currently 
contemplated for building In space becoming 
larger, Identification of the dynamic characteris- 
tics of these structures is an important area of 
research. Accurate on-orbit identification 
becomes a necessity sa these structures cannot be 
assembled fully on the ground becauec of lte size, 
and also It is difficult to predict an accurate 
model on the ground, As the performance require^ 
mencs for these structures In space become 
stringent, however, It becomes imperative to 
Identify their charecteriscice on-orbit and modify 
the control system as dictated by the Identified 
control system design model (adaptive control). 
This paper highlights the model determination 
phase of the adaptive control problem. This phase 
Involves not only determination of parameter 
estimates for sn assumed linear form, but also Che 
order of the linear model form. 

An approach for identifying the dynamics of Large 
Space Structures (LSS) thee estimates model order 
in addition to model parameters Is presented In 
Sundararajan and Montgomery (1983). It uses 
lattice filters which provide an order as well «s 
time recursive algorithm for linear least square 
signal estimation. A comprehensive tutorial on 
Che theory and applications of lattice filters has 
been given by Freldlander (1960). The main 
results from the paper of Sundararajan and 
Montgomery (1983) are! the least square estimate 
of the measurement sequence! the model order 
required to fit the measurements; the associated 
lattice model (this includes mode shape estimates 
that are orthonormal in the measurement space) and 
the associated auto-regressive moving average 
(ARMA) model of the measurement sequence. The 
mode shapes obtained by the lattice, filter are not 
the "natural” modes but a linear combination of 
them. In order to compare the identified mode 
shapes to chose predicted by finite element analy- 
sis, s decoupling method Co obtain natural mode 
shapes from the lattice mode shapes have been 
developed in Sundararajan and Montgomery (1982). 
Using the above spectral decoupling method to 
obtain natural mode shapes, mode frequency and 
damping can be obtained using an equation error 
parameter Identification method (Jol'inaon and 


Montgomery (1979)) chat employs a second order 
ARMA model Co represent the natural mode ampli- 
tudes, This procedure la followed herein. 

The objective of this paper is to present the 
experience In using lattice filcer theory for 
Identification of structural dynamics of two flex- 
ible structures. They consist of a one-dimen- 
sional free-free beam and a two-dimensional 
flexible grid apparatus. The structures are part 
or an experimental facility at the Aerospace 
Control Laboratory at NASA Langley Research Center 
for studying advanced control concepts for large 
space structures. The beam apparatus provides a 
simple structure to test the basic concepts first, 
and the grid apparatus provides a more complex 
structure close to the real spacecraft. Before 
presenting the results of lattice filter identifi- 
cation for these structures, a brief outline of 
the basic theory Is given. Results are presented 
for Che identification of the dynamics of the beau 
using experimental data. Next, the same Is 
repeated for the grid apparatus. Conclusions 
based on the above study results are then 
summarized. 


Summary of Lattice Filter 
Identification Theory 

For the application considered here, we assume 
that the kth measurement sample Is of the form 

y k " [V k) > y 2 (k) y NS (k> ] 

where NS represents the number of sensors. It is 
assumed Chat y Is generated from a model system 
wherein 

y k’% + n k <° 

Here, 9 Is an NSxNM mode shape matrix, q k is the 

NMxI modal amplitude vector, and n|< I s a NSxJ 
gauaslan random variable with zero mean and s 
covariance matrix R. Typically, for structural 
dynamics applications, each component of q^ is 

the output of an uncoupled second order process. 
The Cask here is to estimate the order (NM) and 
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obtain th« least square estimate of q^ fro* tha 
N*I moaeurement sample* yg through y N . 

Sundararajan and Montgomery <199!) praaanc a deri- 
vation of tha equations that re. at* ordar, n, and 
time, I. recuralona for tha normalised forward and 
backward raalduala aa wall aa tha leaat aquara 
estimate of tha itaaauranant vactor. Thcaa equa- 
tion* are Hated below; 

•i,n+l " <l * k i,n+J ) ,/2< ~l,n " k i,n+l£i-l ,n J 

ii,nt| " tl “ k2 l f nai>“ l/2 (. r .l-l,n " k i,n+lil,n> 
N-l 

y H ■ l Q €%, n |i*-r n > < 2) 

wherein 

k l,n+l ’ ^i.n' ^-l ,n* 

and£(x|y) la the othogonal projection operator of 
tha vector x onto tha vector y. The ayametry of 
the recuralon formula* are apparent. Tha equa- 
tion* are coupled by the tern k i((1+ | which i* 

customarily called the "reflection coefficient.” 

The etructure of thi* equation i* depicted in 

Pig. I where w# have used the aynbol r“* to 

repreaent the cine ahift operator, i.e. 

*”* r. - r, . 
i,n 1-1, n 

It should be noted at this point that the lattice 
filter la a modified Gram-Schnldt procedure 
Involving both forward and backward residuals 
wherein the backward residuals form an orthogonal 
basis for the entire observation sequence. Hence, 
any least square eatinate 1' Che orchogonal 
projection onto this basis. Aasunlng at Chi* 
point the order NM has been obtained (which is 
explained below), the lattice filter ha* 
decomposed Che estlnaclon of y into the form of 
equation (1) , 


nodal order between 0 and N+l inclusive. Thi* 
information provide* the basis for the nodel order 
determination me hod described next. 

Clearly, in this approach one nay "fit tha noise" 
by continually increustng the order of the system! 
however, once eh* order of tha estimator has 
Increased beyond the correct order, then the 
residual errors should 11* within * noise band 
which can be predicted based on assumed noise 
characteristics. A threshold value can be 
selected based on this predicted noiae band, and 
order can be determined by a test of whether or 
noe the residual* have been reduced to lie within 
the noise band. Also, the teat can be made 
considering several data sample* when using a data 
window. Sundararajan and Montgomery (1983) docu- 
ment the experience in choosing the data window 
site HU and the threshold level baaed on simula- 
tions. Having deflnrd the order required to fit 
the data using a linear model, for comparison with 
finite element analysis predictions, we seek a 
fixed sec of basis functions that are spectrally 
decoupled. A method to obtain the decoupled modes 
from the lactic* filter modes using digital 
Fourier transform (OFT) has been presented in 
Sundararajan and Montgomery (1982). Essentially, 
at this point we have estimates for order NM, mode 
shapes 0 and modal amplitude time series q(k) from 
the lattice filter, 

Since the ultimate objective of Identification is 
for control system design, an ARMA model is 
Identified using the modal amplitude time series 
q(k). The method is based on an equation error 
method described in Johnson and Montgomery 
(1979). For each mode, the model is described by 
the equation; 

q(k-l) - Aj q(k-2) + A 2 q(k-3) + BjU(k-2) + B 2 u(k-3)’ 

( 3 ) 

The equetlon error is given by; 
e(k-l) - q(k-l) - q(k-l) 

- q(k-t) - (Aj q( k-2) + A 2 q(k-3) 

+ B ( u(k-l) + B 2 u(k-3)} 


wherein 


♦ 


lr 


1,1 * 


Vnm 


I 


and 


■^.NM 


q(i) 


i.e., the backward realduals rj >n form the 

orthonormal basla, or the mode shape raacvlx f and 
the forwerd residual ^ „ repreeente the model 

amplitude time aeriee. 


The lattice filter has the following advantages; 


1. Given e baeia for order N, e basis for 
order N+l can be obtained uelng the 
recuralon formulae. 


2. Becauee of the modified Gran-Schmldt 
procedure, the basis for alL orders n 
between 0 and N are the first n elements 
of the basis of order N. 


3. The estimate assuming any order n 
■oeween 0 and N+l can be computed using 
e. melon (2). 

Tbs lac ice filter provides the Information 

nee. set) lecermine the residual sequence for any 


where q la the model amplitude estimated by the 
lattice filter, u Is Che modal control force, k is 
the (ample number and Aj , A 2 * B [ . B 2 * re 

Che ARMA coefficients. The term in brackets la 
Che model equation. The ARMA coefficients ere 
then updated by; 


A, (kj 


Ajfk-I) 


q(k-2) 

A,(k) 


A,(k-l) 


qCk-3) 


M 


+ o(k-l)p 


Bj(k) 


Bj(k-l) 


q( k-2) 

B,(k) 


B,(k-i) 


q(k-3) 

_ 2 


L 2 


U' J 


The weight p assure* stability if 


0 < U < 2/|q 2 (k-2) + q 2 (k-3) + u Z (k-2) + u 2 (k-3)( 

This identifier performs well in a low noise 
environment, but when the information content of 
the signal la small, it attempts to fit the noise 
(Thau, et. al. (19S2). Also, the ideal ARMA model 
for the beam has Input parameters (B's) which are 
three orders of magnitude smaller chan the (A'a). 
Thle causes a very high sensitivity Co noise In 
Che identification of the B's, and when the Input 
force la applied, it tends to cause Che identifier 
gain on the A's to decrease significantly. 
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Al chi igh these effects are evident in the reiulte 
presented here, they did not prevent successful 
Identification, 

1 f one le interested in determining the deeping 
ratioe and natural frequencies of the node*, they 
can be obtained in a straightforward Banner Iron 
the equation (3), However, it should be noted 
that this is not unique due to the foldover 
phenomenon due to sampling, by finding the roots 
of che equation (3) and uliing the relation 

t * e* T in the primary strip, where t it the 
sampling period, the damping ratio and natural 
frequency can be obtained. 

The bfhevior of this overell system identification 
methodology with experimental data *# discussed in 
the eubaequent sections. 


Experimental Studies for the 
Identification of a Free-free Beam 

In this section, the lattice filter theory 
developed earlier ie llluetrated for the Identifi- 
cation of a one-dimensional free-free beam, The 
identiticacion scheme yields che structural 
dynamic characteristics of the beam. The experi- 
mental epperstui for the free-free beam ie shown 
in Fig. 2, It coneiata of a 12-foot beam of 
rectangular c?o8*-eection which la suspended from 
the cel ting by two cables and ie attached to four 
electromagnetic force actuators. There ere nine 
noncontacting deflection sensors that measure the 
translational deflection of the beam. The actua- 
tors are compensated to eliminate the effects of 
friction as much as possible. This compensation 
is nonlinear, producing a force in die direction 
of the beam motion at Che actuator attachment 
points which is designed to equalise the effect of 
friction. Testing was done manually esdting 
the bean approximately in its first flexible mode 
end sampling the nine sensor s at 66 samples per 
second. A total of } seconds of data was stored 
on a tape which waa poat processed with the 
algorithm, Figure 3 shows a time hietory of some 
of the measurement data processed by the algo- 
rithm, The innovations sequence for sensor A, 
INOV^, is shown Just below its time hietory, 

Also shown is the norm of the forward estimation 
residual, ENORM, which includes all components of 
the measurement vector, Below che norm Is the 
estimate of model order. This was obtained using 
a data window of eight samples. Initially, the 
order escimacor fills the data window, and hence, 
Che Indicated order estimate Increases to 11. 
After this the order estimator settles to 2 
indicating that, even though we attempted to 
excite only one mode, there were, in fact, two 
Significant nodes excited. Note also that the 
norm of the forward estimation error ia email 
compared with the value at the start of che 
process when Che order estimate was settling. 

The modal amplitudes obtained from the lattice 
filter are spectralJy decoupled, using the proce- 
dure discussed earlier, after enough data are 
taken to accurately cake the DFT C 64 time samples, 
about l second). This occurs at about 1,75 

seconds, che first .73 seconds being used for the 
identification of mode shapes and model order (see 
Fig. A). Figure A shows the modal amplitudes for 
both of the identified modes. These are the 
signals that are inputs to the parameter Identifi- 
cation scheme used to identify the parameters of 
the ARMA model of the modes. The Identified ARMA 
parameters are ohown on Fig. A for each of che two 
modes identified. The a priori parameter esti- 
mates are initially offset from the values 
predicted by a finite element analysis which are 
also indicated In Fig. A. These parameters track 


the instantaneous value required to minimize the 
output error. Oru* possible explanation of fh* 
oscillatory behavior of the node 2 parameter 
estimates is the nonlinearity of the actuator 
compensation. Nonlinearity la apparent in the 
aenaor 6 dots on Fig. J, Note that lattice filter 
produce# a linear Isaac aquare fit of eha data to 
the measurement*, and in ao doing, produces a 
predominantly linear first mod* estimate and lump* 
the nonlinear dynamics into tho higher modes. 
Thus, the parameter tracking is more ecable in 
mode 1 and produces estimates of an undamped 
< A2--1 ) oscillation at neaily 2.7 Ha, If the 
algorithm la eonatralncd to an order estimate of 
one, the predominant response ie llnqar, however, 
the fit error le increased by «n order of 

magnitude. 

The mod* shape* estimate* obtained from the 

lattice filter ere shown in Fig. 3, In this 

figure w* compare the estimate* obtainad by three 
method#, one analytic, snd two experimental. The 
analytic result la ch« primary mode ship* of the 
bian using Euier-Bernoulli theory. The two 

experimental result* which are in substantial 
agreement are the nonllnaar laser equere* 
algorithm of Thau, at al. (1982) and the lattice 
filter algorithm of thla paper, igsin note that 
there t* apparently an effect of the four attached 
actuators on (he dynamic* of the test article, 
The lattice filter produces two modes, one near 
the mode of Thau et *1. (1982) and another Chat Is 
shown on Fig, 5. This other estimated mode does 
not resemble any mode ansiyticaliy predicted using 
linear Euier-Bernoulli theory, rather, is required 
to model che effect of nonllnesrlties in che 
apparatus, 


Experimental Studies for the 
Identification of * Flexible Grid 

Next, che lattice filter identification scheme Is 
tested in a more complex structure compared to 
that of the bean. The candidate structure 
considered is that of a two-dimensional riexlble 
grid. Identification results are given using thu 
experimental data obtained from the laboratory 
apparatus. 

Figure ft shows the flexible grid experimental 
apparatus in the Aeroapace Control Laboratory, at 
NASA Langley Research Center. The grid la a 
7 ft x 10 ft planar structure made by overlaying 
aluminum bars of rectangular cross section. The 
bars are centered every foot ao that there aro 
7 vertical and 11 horizontal bars. As shown in 
Fig, 6, the grid is suspended by a cable at two 
locations on the top horizontal bar. The motions 
of the grid perpendicular to the plane of Fig. h 
arc the ones of interest in this study. There are 
nine noncontacting deflection sensors mounted on a 
back frame which give a 9x1 measurement vector. 
The sensor data is linked to the main Cyber 179 
Real Time Computer System at NASA Langley Research 
Center so chat che Identification can be carried 
out In real time. 

In order to compare che experimental results of 
lattice filter identification of che flexible grid 
facility with predicted values, a finite element 
anlayils of the grid was made which Included Che 
suspension cables. Four degrees of freedom appro- 
priate for motion normal to the plane of che grid 
were considered. No damping was Included in the 
model. 

Experiments were conducted uelng the grid 
apparatus described above. The procedure for 
conducting the experiments was to excite the grid 
using an air shaker. The shaker was capable of 
periodically exhausting a Jet of atr that Impinged 
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M I'M |(l< it MUMI IMKIM I K«H* Ml M IM 
M(IM loll hand (•mar «f tha grid. Thd 

tiniiiMf at (la Jat Ml i<|«mtli Iraa 0 U 

M Mi. Tha raawlllng fill radiation una art 

purely IlMIlUll Nl Ml (ltd I* Mraulll. 
M< luN if (M I mil llilll al llM MIImIIM 
••mi an * 0 la appro* iMtaly I la. * (M aaalaua 
paah *• mm dallaatlana al tha grid vara Hatted 
(a Mmi I la. Ida pen* ta paid aallactlaa 

■Mi tail Hall, cha ill ahaaar mi turned all 
art (la (ill waa allowed ta viArata l|N|| Mil 

only air uM aatarlal leaping, a Cykar I/I rail 
llaa caapular lyilaa aenplad tha da I tact laa aaaaar 
lata •’ >2 -Mi. tar l iac . Tha lata vaa it oral an 
« lyilaa lata Ilia tar further aaaiyali. Thu 
tail procedure ••• rapaataM lar aavaral ihahar 
I raaaanc ta tha range al I ta 10 Mi. Tht • 
range haa an upper Halt bacauaa al tha aaapllng 
frequency 02 Ml.) Tha lunar Halt n aalactad ta 
include tha praMIctaM lovaat vibration am 
(2.2 Nr.) light lata lat • carraapandtng ta 
dllfarant ihakar aicllktlaa fraqueaclaa wara 
Mllli an! itaral oa tipaa. Tha following 
dlacuaalon partilna ta raaulta attracted Iroa lata 
••( ll»a. 

Plgure / praianti lata Iroa union I and A ia 
wall aa tha lacttca flltar orlar iitlaata. lat 
thli cm tha order aaciaaia mi kaaaa an a art • 
alnaau al A i«npl«« an4 ipactral licwpling vaa 
aana with k* t laa iiapiaa. Proa tha ftgura, tha 
orlar aatlaata ta iaan ca converge ta an oactlle- 
t Ion kacooaan 2 ana ) at about ,t gMa At kfcoat .• 
tac, tha orAar aatlaata vaa l iaa at 1 ana aata 
col tact mn al tha AA t laa aaapiaa required lar tha 
uPT vaa atartaa at 12 Ml. Tha DfT vaa accoa- 
ptluhod at about 2.1 aac ana than tha Uncoupling 
t ranaloraat Ion antrU vaa calculalaU. Tha aoUat 
aapiltudei altar thta t laa ahaulO contain a alngla 
fraquuncy ana tha tranaloraaU aoUa ahapaa ahoulU 
carraaponU ta tha natural noaaa al tha atructura 
which vara aicItaU. In that aannar, thraa node* 
vara aatractaU Iroa aapartaantai data tapa t. 
Thaaa ha .a traquanciaa naar ,1N|, 2.) Mr., ana 1 
Mi. Tabla I praaanta tha aoUa ahapa aatlaataa 
ohtalnaa Iroa tha aaparlaant. Alao praaantad ara 
•a lac tad aoOa ahapa praUtctlona taken Iroa Itntto 
olvaont onalyooo. Tha ao'oa aalacraU vara thaaa 
vhoaa fraquanclae brackat tha oapartaantally 
durivod mat. Tha lallowtng dlacuaalon daala with 
tha Tabla I data In arUar of tncraaalng 

frequency. A good coapartaon doaa not talk! with 
• tthar brackatlng Malta alaaant analvala aodi ana 
tha flrit aapartaantai aoUa. AUUir lonal ly , thara 
la aoaa banding In cha aapartaantai aoda aa la 
evidenced by tanaark *, I, ana k. Tha aoUal 

veplltud* lar thta ooUi la ah own tn Pig. I along 
with thn AAMA paraaatata A, and Aj lor tha 

aoda ana thatr priaary atrip counarparca al 
Uaaptng C and frequency a. Par tha aaconU aoda, a 
good coapartaon doaa aatar katwan It and tha 2.0/ 
Mi llntta alaaant analyala aoda. Mata that, 

hovavar, tha output ol tha aanaor A la oppoatta tn 
atgn and raducad ta aaplltuda Iroa cha Mntta 
alaaant predlctU.i. Thta aaana that a faadback on 
that aanaor kaaad on tha Mntta alaaant analyala 
util ha daacahlllgtng naar tha 2.2 Ml frequency. 
Tha validity al thta deduction can ha ••tabllohaa 
by anaatnlng tha output* ol aa.iaora I and A 
(Pig. /)• According to cha Malta alaaa.«c analy- 
ala, tha 2.2 Mi contant ol tha aanaori ihould ba 
oppoatta In algn. Mowavor, they ara la phaaa ta 
agruaaant with tha taantlllcatlon raaulta. 


naoaly that al a fraa-fra* bean and a llailkla 
grid atructura. Tha raauita inauata that tha 

lattlca M I lar can ba allacttvaly uaad lar on- 1 1 na 

Identification al tha Makar ol aadaa, a od a 
ahapaa, nodal leaping ana art a I fraquanctat Iroa 
[ha aaaturaaant data. Tha eaperlaentel raaulta 
ala# indicate that thara la caaaldarakla diaagtaa- 
•ant katuaan than and analytical pradlctiona. 
la art an thaaa aapartaantai atodtae. tha aain 
conclualan :hat can ha drawn la that lar largo 
apace atructura*, an-arklt teat mg and IdantlMca- 
t laa ta aaaantlai kalara central I* aitaoptad. 
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Toklo I.- Coaporlaon Satwavn tho Plnlto tloaont 
Prediction* *n3 fa rporTaontoI 
Idontl heat Ion 


Mod* | Nod* 2 Mod* ) 

Caspar laona Cooper leone Caspar loom 

Origin PIPPCPPIP 
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-.2k 

.2* 
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P - Prediction baart on Unit* alanant ar.alyala 
t - Calculation kntod on Idant 1 1 (cat Ion Iroa 

aapartaantai data 


CONCLUSION 

Tha application ol laaat aquar* (attic* filter* In 
ldantllytng tha dynanlc charactarlatlca ol highly 
llartbl* atructura* haa baan praaantad. Tha 
theory haa baan uaad to identify the atructural 
charactarlatlca ol two aapartaantai hardware. 
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ABSTRACT 

Thla papor ravlaua tha uaa of the loaat square 
lattlca filter In adaptive control systems. 
Lattice flltere have been used primarily In apaoch 
and algnal processing, but they have utility In 
adaptive control becauae of thair order-rocuralva 
nature. They are eapectally uaeful In dealing with 
structural dynaelca system* wherein the order of a 
controller required to daap a vibration te variable 
depending on the nueber of node a a 1 gn 1 f leant ly 
excited. Application* are preaented for adaptive 
control of a flexible beam. Alao, difficultly In 
the practical laplementat Ion of the lattlca filter 
in adaptive control are dlacuaaed. 


INTRODUCTION 

For large flexible apacacraft, deaign eodela will 
probably not be adequate. Hence, an adaptive 
control ayatem la highly dealrable. Early reaearch 
into adaptive vibration control of large flexible 
xtructures la reported In reference I. Therein, 
adaptive control of a aplnnlng annular momentum 
control device ( AMCD) waa atudled. That acheiae 
conalated of almultaneoua Identification and 
control with the objective of reguletlng the out- 
of-ptane deflectlona of the aplnnlng AHCh. Some of 
Che dlaadvantage* of the method were the require- 
ment of aelectlng 'he number of mode* to be uaed 
for controller deaign, the uae nf analytically 
predicted mode shapes, and the coupling between 
modes due to lnhomogenltlas In the system. Lattice 
filter adaptive control la a new method which 
attempts to overcome these problems. It la, hence, 
well suited for tiie adaptive control of flexible 
spacecraft. 

The least square lattice filter haa been uaed 
extensively In the field of speech and signal 
procesalr -,2 (reference 2). In these applications 
the filter la designed based on a predetermined 
estimate nf system order. Reference J la a compre- 
hensive tutorial on this subject. Concerning 
adaptive control, reference A proposes a self- 
tuning controller configuration ualng lattice 
filters. Thla scheme requires computing the poly- 
nomial coefficients for the plant and controller at 
each Iteration and enforcing a known feedback 
structure for the controller. Reference 5 propose* 
Inverting the transfer function of the plant for 
general adaptive control. Thla Idea, with the 
least mean squares ( LMS) algorithm, waa utilized In 


reference A to obtain adaptive control. 

Reference 7 proposes a similar approach ualng 
lattice forma Instead of the LMS algorithm. 
Reference A takes this approach but uses a lattice 
model Instead of an autogresslve , moving average 
with exogenous variables (AJUtAX) model where 
familiar controller techniques could be used. All 
of these acheaea attempt simultaneous identifica- 
tion and control or direct adaptive control. For 
each caae stability questions are not resolved 
analytically; neither are simulation results 

available In the open literature. 

As opposed to simultaneous Identification and 
control, the scheme discussed herein consists of 
conducting tests to obtain a design model, vali- 
dating the model, designing a controller based on 
the validated model, and finally, engaging the 
control system. This approach la Ideally suited to 
tha control of large flexible spacecraft because of 
the pesslve environment of outer space and the 
potential of relaxation to a controller that Is 
known to be stable - that of collocated rate feed- 
back. It waa originally presented In reference 9 
and represented the first use of a recursive vari- 
able order structure for adaptive control. 
Therein, the lattice filter was used to provide an 
on-line estimate of the system order, mode shapes, 
and modal maplltude* to provide a validated modal 
control design model. After the Identified model 
parameters are validated through a series of test 
procedures, they are used in a modal pole-placement 
ontrol law design. Figure 1 shows the adaptive 
control scheme using lattice filters. 



Figure 1.- Maptlve Control with Lattice Filters. 


28 


the purpose of this paper Is to assess progress in 
using lattice filters In adaptive control of 
flexible spacecraft and to highlight problem areas 
for further research. First, lattice filter theory 
and order determination is summarized following the 
original development of reference 10. Then, their 
use in adaptive control is dlscursed along with 
applications to the vibration control of a beam,, 
Finally, difficulties arising in the practical 
implementation are discussed. 


suwkart or lattice filter theort 

AND ORDER DETERMINATION 


For application considered herein, we assume that 
the 1th measurement sample is of Che form 

* [y t (0, y 2 <0, .... y NS (D) 

here NS represents the number of sensors. It is 
assumed that y la generated from a model systeu 
wherein 

y i ” ** T t + v i <l> 

Here, t is an NS x NH mode shape matrix, ^ is 
the NM x 1 modal amplitude vector, v^ is a 
Gauss lan-random variable with zero mean and 
covariance matrix R. NM represents the number of 
modes in the system or order of the system. 


Reference 10 presents a derivation of the equations 
that relate any order, n, and time, i, recursions 
for the normalized forward and backward residuals 
as well as the least square estimate of the 
measurement vector. These equations are listed 
below: 


— i ,n+l 
-i,n+l 


(1-k 
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i,n+l 


(1-k 
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( -i,n" k i,n+l-i-l,n } 

( ~i-l,n _k i,n+l-i,n* 


wherein 

k. “ <e. , r > 

i,n+l — i,n — i-l,n 


and < > represents an inner product. The symmetry 
of the recursion formulae is apparent. The equa- 
tions are coupled by the term k£ >n+ j which is 
customarily called the "reflection coefficient." 
The estimate of the measurement (reference 10) at 
sample i for a model of order n is 


*».»■ 'h/ 

where (o represents an orthogonal projection 
operator. Hence, 


y i,n " l r n-1 ,0 ,r n-l ,1 1 


n-1 ,n-l 1 


so that 


1.1 

V 

L k i» n J 


y i " V T i + ^i,n 


( 2 ) 


where is an orthonormal NS x n matrix 

(ln-1,0 In-1, n-l J generated from the 

lattice filter, and la the n dimensional 

vector of reflection coefficients and jjj n is the 
NS dimensional estimation error vector. ' 


Clearly, in this approach one may "fit the noise" 
by continually increasing tha order of the system; 
however, once the order of the estimator has 
increased beyond the correct order, then the 
residual errors should lie within a noise band 
which can be predicted a, priori based on assumed 
noine characteristics. *~A threshold value can be 
selected based on this predicted noise band and 
order determined by a test of whether or not the 
residuals have been reduced to lie within the noise 
band. The residuals will generally consist of 
signal and noise parts - the signal part being 
reduced as the correct order is reached until the 
residuals essentially consist only of noise. This 
test is carried out based on a data window of NW 
samples. Thus, assuming that the data can fit a 
linear model and chat Che noise process is 
Gaussian, for i large enough, 


NW T T 

E( l ej „ e. ) - NW E(v v!) 

tii l » n ‘.h 1 1 


t ; 1 

- NW tr E(v.v‘) - NW \ o~ (3) 

j-1 J 

where E is the expectation operator. This can be 
used as the one sigma threshold for the order 
determination test. In the last equation oj is 
the standard deviation of the noise process 
for the Jth sensor. Reference 10 documents experi- 
ence in choosing the data window size NW and the 
threshold level based on simulations. 


ADAPTIVE CONTROL USING LATTICE FILTERS 

Independent Modal Space Control ( IMSC) 
(reference U) is a control scheme specifically 
designed to deal with flexible spacecraft in a 
modal form amenable to control law design. 
Unfortunately, it requires natural modes and not 
the oithonorraal basis provided by the lattice 
filter. Consequently, in order Co interface the 
lattice filter outputs with Che target adaptive 
control scheme (figure 1) and to make comparisons 
with finite element analysis predictions of natural 
modes, a method is needed to obtain natural mode 
shape estimates from the lattice filter basis. The 
filter updates tha NM basis vectors at every sampls 
instant. While the order estimate NM remains 
constant, the updated basis vectors are related by 
a mere rotational transformation. The assumption 
of che target adaptive control scheme is that the 
system motions can be modelled by a constant and 
finite set of natural modes and their associated 
modal amplitudes over a reasonably long time 
interval. Hence, when the estimated system order 
is constant, the basis elements used to derive the 


model amplitude time aeries required by the target 
adaptive control scheme are not changed* 


e(i-l) - r^(i-l) - [Ajd-Dyi-2) + A 2 (l-l)T N (l-3) 


The transformation from the lattice filter to a 
natural mode basis should satisfy 


+ Bj(l-l)f(l-2) ♦ B 2 (l-l)f(l-3)] (6) 


wherein the subscript L refers to the lattice 
filter and N refers to the natural modes. A non- 
singular matrix T, will satisfy this condition 
provided 

V ■ V» ■ w 

Since the order estimate Is assumed constant, this 
matrix can be approximately determined on-line 
using the digital Fourier transform (DFT). Herein, 
this Is accomplished aB follows, Since the lattice 
fitter uses the current measurement sample as Its 
first basis element, the corresponding modal ampli- 
tude time series contains NM frequencies. Hence, 
the DFT spectrum of this series will contain NM 
peaks corresponding to these frequencies. The 
frequencies (wj, i» 2 , ..., can thus be 

Identified by searching this spectrum for these 
peaks. Assuming that the motion la comprised of 
undamped structural vibrations, the matrix T, which 
produces the desired transformation can be calcu- 
lated as 


*•(*[<«!>] 


"K'Vl 


wherein, [^(w), ^ (“)) io a " NM dimen- 

sional vector of the modal amplitude transform. 
Using thiB matrix, the digital Fourier transform of 
each component of Will be zero at the discrete 
frequencies, w . , jvi. :)»* item which degrades 
this approximation is the error in using DFT 
instead of the true Fourier transform. Still 
another Is the assumption that the motion is made 
up of undamped structural oscillations. In spite 
of these items, reference 12 shows that this 
approach produces good estimates of the natural 
modes for the beam used herein. 


The decoupled modal amplitude time series, 
as obtained above in equation (4), is then ana- 
lyzed, for each mode, to identify the parameters of 
its autoregressive, moving average (ARMA) model. 
The Inputs to each ARMA modal model are the 
generalized forces and hence, each model takes on 
the form: 


The gradient technique of reference 1 is used to 
identify the parameters p ■ (A*» B^) using the 
iteration sequence 

p(i) - p(i-l) + e(i-l) [WjT N (i-2), V 2 yi-3), 

W 3 f(i-2), W 4 f(i-3)] (7) 

As indicated in reference 1, the welghta W must be 
selected consistent with the relation 

HjT^l-2) + W 2 tJ(i-3) + W 3 f 2 Ci-2) 

+ W 4 f 2 (l-3) < 2 

and the inputs to the algorithm, and f, must 
be sufficiently varying and large If the parameters 
are to converge to their correct value. 

For the identification and control scheme explained 
above to work satisfactorily in a closed loop 
environment, it is necessary to validate the design 
model. Three tests are suggested herein which 
check the following: 1) model fit error; 2) 

parameter convergence; and, 3) signal informa- 
tion. These tests have been used successfully in 
simulation and experimental work. The fit error 
test uses a fixed parameter set to calculate an 
estimated modal displacement for the past NT 
samples. 


°fit > £ V*" 0 ) " {Aj^d-n-l) + ^^(1-0-2) 
n**0 

+ Bjf^i-n-l) + B 2 f(i-n-2)), k > NT 

If the absolute sura of the error between the 
modelled displacement and the displacement calcu- 
lated by the lattice filter exceeds a given thresh- 
old, the fixed parameter set is updated with the 
present identified parameter set. This process Is 
repeated until the parameter set fits the data. 
The convergence test runs concurrently with the fit 
test. It simply checks the magnitude of the 
changes in successive estimated parameters. 

NT 

O > I |p ” P I for P T “ [a, , A, , B, , B,J 
qonv a I ‘n n-1 1 1 1 2’ 1’ 2 J 


V‘> ■ -SV 1 -” * W‘- 2> 

+ B + B fCl-2) (5) 

where the f represents the generalized forces. 
Given the and f's, the parameters A and B 

above are Identified and used in the control law 
design process. Thus, the ARMA model output error 
is 


If the absolute sum of ten successive parameter 
estimates changes is above a specified level, a 
logical switch is set to indicate failure. The 

third and final test is on information content of 
the estimated modal amplitude signals from the 
lattice filters. The purpose of this test is to 

check whether enough information is present in the 
signal for proper Identification of the 

parameters. If this test falls, the controller 
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gains* are not updated based on the identified 
parameters, but are frozen at the last values 
before the teat failed. Here, the estimated modal 
amplitudes and velocities from the lattice filter 
are checked for sufficient excitation by summing 
over ten samples. 



The second term in the above equation represents a 
measure of velocity estimates. If the sin is belov 
a threshold, o <nf , the updating of the control 
gains based on the identified parameters is 
stopped. The information and fit error tests 
constitute one test for each mode and the conver- 
gence and reasonability teats constitute four tests 
for each mode. Thus, six tests must be passed 
before control is applied to a given mode. The 
actual stability and performance of the controller 
is directly affected by the criteria chosen for 
passing a test. If the test criteria are too 
stringent, system noise and noniinearlties may 
preclude initiation of control. However, If the 
teats are not adequate, it la possible that an 
error in the estimated parameters could result in 
gain calculations which produce an unstable 

system. 

Now, consider the philosophy to be used when the 
tests described above pass or fall. When alL the 

tests for parameters of a given mode have passed, 
control gains are calculated according to a 
previously developed pole placement scheme 
(reference 1). The control force commands are then 
calculated using these control gains. Considering 
the philosophy used when the tests fail, two cases 
were studied. In the first case, chen the tests 
failed, control was turned off and the control 
forces were made zero. In the second case, when 
the tests failed, updating of the control gains was 
stopped and they were frozen at their values prior 
to the test failure. In this case, the control 

forces were not made zero and were computed using 

the frozen control gains. From a detailed study of 
both cases t It was found that the performance of 
the adaptive control system in Che first case was 
superior to that of the second case. 


APPLICATION TO A FLEXIBLE BEAM 

The closed-loop adaptive control scheme of figure 1 
has been tested in the digital simulation for the 
12-foot, flexible free-free beam located at NASA 
Langley Research Center. The simulation contains 
the mathematical model Of the beam apparatus in 
modal form. For this study, the simulation 
contains one rigid-body mode, Che first three 
flexible modes, nine deflection sensors, and four 
actuators for control purposes. The initial 
conditions on Che modal displacements were 3et to 
.05 in. and the modal velocities were set to zero. 
The modal damping was also set to zero. A digital 
sampling rate of 32 Hz was selected for the simula- 
tion, and the standard deviation for all measure- 
ment noise was assumed to be .005 based on observed 
noise In the available hardware. The lattice 
filter estimates were based ot: a data window size 4 


(reference 10). The testing procedures were all 
carried out based on data window (NT) of ten 
samples. Initial parameters estimates were offset 
from the mathematically correct values to test and 
verify the rapid convergence of the identification 
algorithm. An arbitrary delay of 2 seconds was 
added between the time Identification starts and 
when the control would be applied to show the 
behavior of the identification scheme. 

At the start of the simulation, the lattice filter 
determines the number of modes in the simulation 
along with the mode shapes. Modal amplitude time 
hiatorles are then generated, From the lattice 
filter mode ahapeB and modal amplitudes, natural 
modes and modal amplitudes are obtained through a 
linear transformation explained in Che earlier 
section. The application of the transformation is 
delayed for 2 seconds because the online 
transformation technique of reference 12 requires 
2 seconds of data for a digital Fourier transform 
data base to obtain the required transformation. 
The natural modal amplitudes are input to the 
equation-error parameter identifier which 
identifies the ARMA parameters. The identification 
results are chen tested using the test procedures 
described above. When the teats are passed, the 
control is turned on. Results of the simulations 
are presented in figures 2-4. 

Figure 2 shows the estimated modal displacement for 
Che first lattice filter mode. The order estimate 
plot shows that the correct order of 4 la obtained 
in .3 seconds. After the parameter identification, 
when all the tests are passed, the control is 
turned on at 5,5 seconds and the modes are damped. 
The result of the adaptive control on the natural 
modes is shown in figure 3. It is evident that 
when the identification is validated by passing Che 
tests and control turned on, the vibration 
suppression is achieved. When Che modes are damped 
out the lattice filter order estimate drops from 4 
to l indicating the flexible modes are damped out. 
Although the lattice filter order decreased, the 
control design order was maintained at 4 throughout 
the time interval when control was on. Allowing 
Che order to vary in real time and updating the 
control order Is a topic for further studies. 
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Figure 2.- Typical time history of an adaptive 
control run using identification, 
testing, and control design. 
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Figure Tine hlatorlea of thro* natural nodee 
with th* lactic* fllt*r ordar Indicated. 


Th* main result* of th* Identification and th* t**t 
procedure* are summer lied In figure 4. For th* 

flret flexible node, th* fig re shows th* 1 1*# 
hi at or lea of Identified frequency parameter A ( 
th* fit error, a parameter that lndlcatea algorithm 


convergence, and a parameter that Indicate* 
Information content of the measurement a, When all 
the teete are paaaed, th* correepondtng paaa 
parameter (plotted aa a binary logical varalbl*) la 
aat to on*. Th* varloua thresholds for th* teata 
are alao narked to Indicate when th* teata paaa. 
Thee* three ho Ida war* determined baaed on detailed 
senaltlvlty atudle* of the nodal control achene for 
th* bean (reference 11). An error wa* lntentlonlly 
put on th* Initial eatlnat* of Aj eo that th* 
convergence of th* eatlnat** to th* correct value 
could be obaervod. When the Identifier la turned 
on, th* eatlnat* converge* to th* true value of 1.1 
f ran 3. The threehold* Indicate that th* fit error 
teat la paaaed flrat and then th* convergence 
teat. With enough algnal In the neaaurenenta th* 
Information teat la alwaye paaaed. When all th* 

teata are paaaed at 1.1 aec, the control la turned 
on. When control la fully effactlva, chat la Wien 
th* node* are danped out, the neaaurenent data will 
contain only th* nola* and th* Information teat 
will fall. Thla la immediately aeen from th* 
hlatory of Aj ** It atart* oaclllatlng with large 
aoplltud* Indication that th* nodal amplitude 
algnal contain* ealnly nola*. Alao, If th* 
paraoeter excursions are large, the convergence 
teata util alao fall Indicating a failure for the 
binary variable paaa. Once thla happen*, th* 
control gain updating la atopped, and control 
force* were made :,aro. 
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Figure 4.- Tine hlatorlea of the teat variable* for 
one node with the teat threaholda and 
logic aum of th* teata indicated. 


Th* adaptive control achene of figure I t* good 
from th* engineering point of view alnce only vali- 
dated model* are uaed for control ayaten dealgn. A 
natural queatlon arlae* aa to th* coura* of action 
when validation teat* fall. The operating environ- 
ment for large flexible apacecraft la, fortunately, 
benign and a ayaten dealgned to aupprea* vibration* 
can be ehut down at th* expenae of having to 
conduct relatively lor.g tern naneuver*. Another 
aavlng feature of large flexible apacecraft la that 
collocated rat* feedback la atable and relaxation 
of th* ayaten to thla node of operation la alao 
poaelble, again, with corresponding degradation in 
performance. Therefore, two option* that can be 
evoked are; one, to ehut down the control ayaten 
and th* other, to revert to a robuat control ayaten 
dealgn which lnaure* atablllty. 

At flrat glance one nay wlah to uae the ARMA nod el 
generated by the lattice filter directly In the 
dealgn proceas rather than ualng 1NSC with It* 

requirement of generating natural node*. Unfortu- 
nately , the current online dealgn capability for 

controller* of vector ARMA processes la not ade- 
quate. Having aelccted IMSC, one nuat obtain 
natural node* from the vector ARMA model or from 
the meaeuremenr tine aerie*. Here the a a roe problem 
arise*, that la, th* current capability of 

eigenvalue/vector analysis for vector ARMA 
processes la inadequate for online Implementation. 
Hence, a tine carles analysis ualng a DFT has been 
selected. The accuracy of the process of 
extracting natural nodes la directly affected by 
the number of data points processed. Hence, there 
la a tradeoff to be made between the higher 
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complexity in control computation# versus the error 
in the natural modes using the OFT approach# Also, 
significant computational saving results if the 
approximation of zero damping can be evoked, If 
this approximation cannot be made, Chen one must 
work with complex modes. 

Since several approximations are required by the 
system, a method of validating the models used in 
the online controller design is essential. 
Analytic methods of validating models based on 
statistical error analysis (e.g. Cramer-Rao bounds) 
are not adequate. Currently, teats on fit erroc, 
algorithm convergence, information content of the 
measurements, and reasonability have been u*ed. 
The thresholds and design constants for these testa 
can be determined only by exhaustive simulation 
and/or hardware tests and is not an online 
procedure. 

CONCLUSION 

This paper reviews the use of the least square 

lattice filter in adaptive control systems. 
Emphasis is placed on the Integration of the 
lattice filter into a practical parameter adaptive 
control system. One novel feature of the 
recommended system is the inclusion of a design 

model validation scheme based on model fit error, 
algorithm convergence, and signal information 
content. An application is presented for adaptive 

control of a flexible beam. These results indicate 
that the lattice filter adaptive scheme Is 

practical for vibration control of large flexible 
spacecraft. Difficulties in the practical 
implementation of Che lattice filter in adaptive 
control are also discussed. These centered around 
the computational burden of transforming lattice 
filter modes into natural modes and the selection 
of the thresholds for online validation tests. 
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APPENDIX D 


ROBUST CONTROLLER SYNTHESIS FOR A LARGE FLEXIBLE SPACE ANTENNA 


ROBUST CONTROLLER SYNTHESIS FOR A LARGE 
FLEXIBLE SPACE ANTENNA 

k 

N. Sundararajan ,» S. M. Joshi, and E. S. Armstrong 
NASA-Langley Research Center 
Hampton, Virginia 23665 

SUMMARY 

This paper investigates the application of the linear-quadratic- 
Gaussian (LQG)/ loop transfer recovery (LTR) method to the problem of 
synthesizing a fine-pointing control system for a large flexible space 
antenna. The LQG/LTR approach of synthesizing a multivariable controller 
in the frequency domain is selected because large flexible structures can 
be modelled with elastic mode transfer functions as additive perturbations 
on the rigid body model and the LQG/LTR approach uses this formulation 
naturally for robust control design. The study is based on a finite 
element model of the 122 m Hoop/Column antenna, which consists of three 
rigid-body rotational modes and the first ten elastic modes. A robust 
compensator design for achieving the required pointing performance In the 
presence of modeling uncertainties is obtained using the LQG/LTR method. 
For the Hoop/ Column antenna, a satisfactory controller design meeting the 
desired bandwith of .1 rad/sec and ensuring stability with unmodelled high 
frequency modes was obtained using only a colocated pair of 3-axis 
attitude sensors and torque actuators. This study also indicates that to 
achieve the desired performance bandwidth of 0.1 rad/sec. and to ensure 
stability against higher frequency elastic modes, the design model should 
include the first three flexible modes together with the rigid body modes. 


* Old Dominion University Research Foundation, Norfolk, VA 
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INTRODUCTION 


» 

One of the planned activities of the NASA's Space Transportation 
System is the placement in earth orbit of a variety of large space 
antennas. Potential space missions will require antennas and structures 
ranging from 30m to 20km in Bize. Applications include communications 

(mobile), remote sensing (soil moisture, salinity, etc.), deep space 
network (orbital relays), astronomy (x-ray, observatory, optical array, 
radio telescope, very long baseline interferometry, etc.), energy and 
space platforms. Specific missions have been pinpointed and future 
requirements have been identified for large space antennas for 

communications, earth sensing and radio astronomy [lj. Particular emphasis 
is placed on mesh-deployable antennas in the 50-120 meter diameter 
category. One such antenna is the Maypole (Hoop/Column) antenna, shown 

schematically in Figure 1, basically consisting of a deployable central 
mast attached to a deployable hoop by cables held in tension [2J. The 
deployable mast consists of a number of telescoping sections, and the hoop 
consists of 48 rigid segments. The reflective mesh, which is made of knit 
gold-plated molybdenum wire, is attached to the hoop by graphite fibers. 
The mesh is shaped using a network of stringers and ties to form the radio 
frequency (RF) reflective surface. In order to achieve required RF 

performance, the antenna must be controlled to specified precision in 
attitude and shape. For example, for missions such as land mobile 

satellite system (LMSS), which is a communication concept for providing 
mobile telephone service to users in the continental United States, It is 
necessary to achieve a pointing accuracy of 0.03 degree RMS (root mean 

square) and a surface accuracy of 6 mm RMS. It is also necessary to 
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have stringent control on the motion of the feed (located near one end of 
the mast) relative to the mesh. Because of its large size and relatively 
light weight, the antenna is highly flexible, with a large number of 
significant elastic modes. Its dynamics can be represented by partial 
differential equations, or by very large systems of ordinary differential 
equations. The resulting equations have many resonant frequencies, some 
of which may be very low, and possibly closely spaced. The natural 
damping is usually very small. For these reasons, control of large space 
structures is a challenging task (3J. Since the system is inherently of 
high order, a practical controller has to be based on a reduced-order 
"design" model. Furthermore, the parameters (i.e., frequencies, mode 
shapes, and damping ratios) of the’ system are known imprecisely. This 
introduces additional modeling errors. Typically, the modeling errors for 
finite element models increase substantially with increasing modal 
frequency. 

Reduced-order control synthesis for the Hoop/Column antenna using the 
standard LQG theory was investigated in [4,5]. The standard LQG procedure 
yielded satisfactory control, i.e., rigid-body bandwidth of up to 0.25 
rad/sec, satisfactory time constants for the elastic modes, and acceptable 
root mean square (RMS) pointing errors in the presence of sensor noise. 
It should be noted that the LQG approach in reference 4 used a large 
number of actuators and sensors (four 3-axis torque actuators and four 
3-axis attitude and rate sensors). It was found in [4] that the first 
three flexible modes had to be included in the "design" model (in addition 
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to the three rigid nodes) to obtain satisfactory performance. The main 
problem with the LQG method was that a large number of weighting 
parameters had to be simultaneously adjusted to obtain a good design. In 
addition, the stability robustness property with respect to Inaccuracies 
in the modal parameters could not be properly evaluated because it was 
difficult to affectively characterize the bounds cn modeling errors in a 
time-domair. setting. In order to reduce these difficulties, normally one 
checks the control design for robustness after the control design is 
completed using LQG or any other method. Such an approach in the 
frequency domain using singular value measures was presented in [6] for a 
large space structure* using different control design methods like LQG, 


//rormnnmr chonar)* t at n - [J rV 1 ^ V A f Vtf> afjQUA mgthodf? . 
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the LQG/LTR approach provides a means of including the robustness-to- 
uncertainities, in the control design process Itself. Since it is in the 
frequency domain, it extends the basic frequency domain design guidelines 
like bandwith, cross-over frequency, etc. from a scalar system to a 
multivariable system. 

The newly emerging robust control synthesis methodology which uses 
frequency domain matrix norm bounds (i.e., singular values) has received 
considerable attention in the recent literature [7-9]. The basic 
framework for frequency domain synthesis using the LQG/LTR methodology was 
developed in [7-9]. It has been applied to diverse tiystems such as power 
systems [10] and aircraft engine control [11]. The LQG/LTR design 
philosophy uses a low-frequency "design model" of the plant and a 
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high-frequency characterization of the modeling errors. This method, 
which characterizes unstructured uncertainty with singular value bounds, 
appears to be particularly well suited for the control of large flexible 
spacecraft due to the considerable uncertainty that inherently exists in 
the mathematical models. 

The purpose of this paper is to investigate the use of LQG/LTR 
multivariable frequency domain methodology in the design of an attitude 
control system for the Hoop/Column antenna, A low order compensator is 
obtained by treating a sequence of finite element design models ordered 
with increasing modal frequency and chooaing the final design model as the 
first one which allows the perforaance/robustness objectives to be met. 
In this sequence of design models, the first one consists of the rigid 
body modes only. Subsequent design models are obtained by the successive 
addition of flexible modes. The designs use 3-axis torque actuators, 
colocated attitude sensors, and attitude feedback. 

The organization of this paper is as follows: The mathematical model 
of the system is described in section 2. The control objective is briefly 
discussed in section 3, followed by a brief description of the LQG/LTR 
technique in section 4. The reduced order (low frequency) design model 
and the high frequency model uncertainty barrier are also discussed in 
this section. Section 5 presents the results of synthesizing the 
controller based on the above procedure using only attitude feedback. 
Some of the problems and limitations observed are also highlighted. Based 
on the study results, the conclusions are summarized in section 6. 
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2. MATHEMATICAL MODEL 

■ 

As a consequence of its large size and light weight, the Hoop/Column 
antenna is a highly flexible system having a large number of significant 
structural modes. A finite element model of the antenna [Ref. 2) is used 
in this paper. The mathematical model considered consists of rotational 
rigid-body dynamics (about the three axes) and the elastic motion. We 
assume that the control will be accomplished by using n T three-axis 
torque actuators. The linearized equations of motion are: 


I a 
s 


r T 

l Tj 

j-l J 


q + Dq + Aq" 4> u 


( 1 ) 

( 2 ) 


where I s is the 3x3 inertia matrix, Tj is the 3-axis torque applied 

by the jth actuator, a s “C 4> s » 0 R , <j/ s ) T denotes the rigid-body 

attitude, q is the Hq x 1 modal amplitude vector (for nq structural 

modes), D » 2 diag(pj(op P2 bi 2> •••» Pn^n^ Is inherent 

damping matrix, (where is the damping ratio for the ith mode). $ is 

the m x nq "mode-slope" matrix (where m * 3n>p) , u *= (T^, •••» 

T t 2 

T n ,j,) r is the m x 1 vector of actuator torques, and A *= diag(o)|» 

2 2 

ai 2 > •••» Wjj ) where is the frequency of the ith elastic mode. 

The rigid-body parameters and the first ten elastic frequencies are given 

in Table 1. Each value of p^ is assumed to be 0.01 for i— 1 ,2 n . 


40 


6 


* 

Normally, the sensors used Include attitude and rate sensors. A 
3-axis attitude y a at a sensor (e.g. a star tracker) output is given 

by: 

y - a + i|)q + w (3) 

where t|/ is the 3 x n q mode-slope matrix at the sensor location, and w 
is the sensor noise. If an attitude rate sensor (e.g. a rate gyro) is 
used, the sensor output y r is given by an equation similar to (3), 
except that a 8 and q are replaced by a s and q, respectively. Torque 
actuators and attitude sensors are assumed to be located near the top of 
the mast at the antenna feed (Fig. 1.) 

Defining x * (a s , a s , q T , q T ) T an n x 1 vector, the state 
space model can be written in the form: 


x u Aj, x + B F u 

(4) 

y * C_, x + w 

r 

(5) 


The sensor noise w is not used in the design process in this paper; 
however, it will have to be included when computing the RMS pointing 
errors. Ignoring the noise, the transfer matrix between the input (3-axis 
torque) and the output (3-axis attitude) is given by: 

G(s) * Gj(s) 4- &2( s) 


( 6 ) 
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where 

G^b) - I^/b 2 (7) 

n^ 

G 2 (s) « l ( V J)/<s 2 + 2p i w i s + <8) 

i-1 

and ^ represent the mode-slope matrices at the sensor and 
actuator locations corresponding to the ith mode). 

3. DESIGN OBJECTIVES 

The basic design objectives for the control systems are: (1) To 

obtain sufficiently high bandwidth (i.e. closed loop frequencies 
corresponding to the rigid body modes) and satisfactory closed loop 
damping ratios for the rigid body and structural modes; and (2) To obtain 
satisfactory RMS pointing errors, feed motion errors and surface errors. 
The first design objective arises from the need to obtain sufficiently 
fast error delay when a step disturbance such as sudden thermal distortion 
caused by entering or leaving Earth's shadow) occurs. The second design 


objective 

arises from the RF 

performance requirements. 

These 

two 

objectives 

may not necessarily 

be 

compatible, 

and may 

even 

be 

conflicting 

» For example, the 

use 

of increased 

feedback 

gains 

for 


obtaining higher bandwith and damping ratios will, in general, result in 
higher r.m.s errors (because of the amplified effect of sensor noise) 
beyond a certain point. Therefore, it is necessary to carefully consider 
the trade-offs between the speed of response and lower RMS error. In this 
study, the main control system specification is that a minimum bandwith of 
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0,1 rad/sec for the closed loop system is to be ensured. The upper limit 
on the low frequency gain is not specified, but it is desired that it 
should be as high as possible. Also, for this cuudy no specification on 
RMS errors was made and this aspect along with measurement noise will be 
considered in the future. 


h. THE DESIGN PROCEDURE 

The LQG/LTR method has been described in detail in [7-9). Here, the 
main steps are summarized first and then each step is discussed in detail. 

(1) Define a "design” model of the nominal plant which is an acceptable 
low frequency representation. Define the high frequency uncertainty 
(robustness) barrier, and the low frequency performance barrier. 

(2) Design a full state feedback compensator based on the steady state 
Kalman-Bucy filter (KBF) > This assumes that the ‘ loop is broken at 
the output. Adjust the weighting matrices in the KBF design until 
it; ; > frequency response meets the robustness specifications at high 
frequencies and bandwidth specification at low frequencies. 

(3) Design a LQ regulator to asymptotically "recover" the frequency 
response obtained in step 2. 

(4) Verify stability, robustness, and performance for the entire 
closed-loop system. 

The first step, which consists of the definition of the plant and the 
uncertainty (robustness) barrier, is often the most important one. The 
basic problem in controlling a flexible structure is the presence of a 
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large number of lightly damped structural modes. Practical limitations 
necessitate the use of reduced-order controllers. Therefore, the 

uncontrolled modes , as well as the error in the knowledge of the 
controlled modes, represent uncertainty. Since the number of structural 
modes is usually large and finite element modeling accuracy typically 
decreases with increasing model frequency, the design model Bhould consist 
of the rigid-body plus the first few elastic nodes. The remaining 
structural modes then (partly) constitute the plant uncertainty. In order 
to obtain an acceptable low-frequency representation, the design model 
must include at least the three rigid body modes. The uncertainty 

barrier is a measure of the plant uncertainty at high frequencies. The 
plant uncertainty can be represented as either multiplicative or additive 
uncertainty (Fig. 2). Additive uncertainties are of the form 

G' « G + AG 

while multiplicative uncertainties are of the form 

G* - (I+A)G 

Multiplicative uncertainty form is the preferred form in the literature on 
robustness studies as the compensated transfer function has the same 
uncertainty representation as the raw model. However, since flexible 
structure models exhibit naturally the additive uncertainty form of the 
transfer function matrix, this will be used in the following studies. The 
LQG/LTR approach requires the characterization of the uncertainty in terms 
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of a frequency-dependent upper bound. Frequency domain sufficient 
conditions are used to test the robustness in the presence of 
uncertainties within that bound. 

For the case of multiplicative uncertainty Lp(s) of figure 2a, the 
closed-loop system is stable if 

o[l. p (jw)-l] < £ [l + (G p (jw)G c (jm))” 1 ] (9) 

where G p (s) and G c (s) are the design model (plant) and compensator 
transfer matrices, and o and £ denote the largest and the smallest 
singular values of the argument matrix, respectively. At high 

frequencies, assuming R [ L p < Jo>) ] R » 1 and I [G p ( jw)G c ( jcu) ]j « 1, (9) 

approximately yields 

ff (G O ) < — ! (10) 

p « <y 

The "uncertainty (or robustness) barrier" is an upper bound l m ( u> ) on 
o(L p )* The system is stable in the presence of such unstructured 
uncertainties if a[G p G c ] < (to) at high frequencies. 

When the additive uncertainty formulation (Fig. 2b) is used, a 
sufficient condition for stability robustness is given by [12] 


a (I + G G ) 
= P c 


> a (AG) 


( 11 ) 


At high frequencies, assuming IG p G c l « 1, (11) (approximately) yields 

a <G c ) < 1 /a (AG) (12) 


That is, the compensator must roll off sufficiently rapidly at high 
frequencies. The main objective of the LQG/LTR approach is to first 
design a full state compensator (based on KBF) which has the behavior of 
the desired loop transfer matrix (i.e,, the loop gain G p G c ). 
Therefore, (from step 2) any loop shaping should involve the product 
G p G c rather than G c alone as in (11) and (12). Assuming that G p 

is a square matrix, 


or 


G - G" 1 (G G ) 
c p pc 


a (G ) < a (G" 1 ) a (G G) 
c p pc 


a (G ) < a (G" 1 ) a (G G ) 

c p p c 


(13) 


(14) 


Using (12) and (11), the following sufficient condition for stability 
robustness is obtained: 


a (I + G G ) a (G ) 

2-5 2- > a ( AG) (15) 

5 < G p G c> 

The second step in the design procedure is to design a full state 


feedback compensator having desirable singular value properties 


The 


L 
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performance of the closed-loop system depends on the low frequency gain 

and the crossover frequency of the loop transfer matrix G p G c ; that is, 

on the behavior of a [GpG c ]. Larger low frequency gain and crossover 

frequency indicates better tracking performance. Thus, £[GpG c J should 

lie above the performance specification as shown in Fig. 3a. The other 

requirement is the stability robustness in the presence of model 

uncertainties. If the multiplicative uncertainty formulation is used, 

according to (10), the a [ G p G c ] plot should pass under the robustness 

barrier a ^(Lp) at high frequencies (Fig. 3a). On the other hand, if 

the additive formulation is used, the robustness condition (15) should be 

satisfied (Fig. 3b). The advantage of an LQG-based full state design is 

•» 

that it has excellent classical properties, and its frequency response can 
be shaped in the desired manner by varying the weighting matrices [8]. 
As discussed in [7], this design can be accomplished using the LQR 
Riccati equation if the loop is broken at the plant input, or the KBF 

Riccati equation if it is broken at the point where the residual signal 

enters the KBF. Herein we select the latter because the objective is to 

control the attitude output. This selection is also consistent with 
[9-11]. The KBF equations are: 

AE + EA T + LL T - ~ EC C T E « 0 (16) 

H * - EC T (17) 

y 
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where L and y are the design parameters, L being an m Q matrix, and y a 
scalar. The matrix H is the KBF gain and l is the corresponding Riccati 
matrix. The KBF loop transfer matrix is given by: 

Gj^Cs) - C (si - A) ~’ 1 H (18) 

Generally, the frequency response cr(Gjyr(ju))) would shift higher as y 
decreases, and the crossover frequency can be adjusted by changing L [6]. 

Having obtained satisfactory singular value behavior of KBF, the next 
step is to design a LQR to "recover" the desired frequency response. This 
is accomplished by solving the algebraic Riccati -equation 

A T P + PA T - PBB T P +qC T C-0 . (19) 

where P is the Riccati matrix and q is a positive scalar. The control 
gain matrix G is given by 

-1 T 
G - R BP 

It has been proven In references 7 and 8 that the loop transfer matrix 
G n G for the overall system (consisting of the plant, the KBF and the 

pc 

LQR) tends to Gj(p(s) as provided that the open-loop plant has no 

transmission zeros in the right half plane. The compensator G c (s) after 
recovery is given by: 

G c (s) - G (si - A + BG + HC)""*H 
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Since the compensation obtained has no guaranteed robustness 

properties, the last step will consist of testing the eigenvalues of the 

entire closed-loop system to ensure stability and robustness. If 

instability is discovered, it will be necessary to return to step 2 and 

redesign the KBF for lower bandwidth and the LQR for robustness recovery. 

If thiB does not produce satisfactory results, it would then be necessary 

to return to step 1 and include more elastic modes in the design model. 

Application of the foregoing LQG/LTR procedure for the Hoop/Column antenna 
# 

is described in the following section. 

S. CONTROLLER DESIGN BY LQG/LTR METHOD USING ATTITUDE FEEDBACK 

The foregoing procedure has been "applied to the Hoop/Column antenna 
model. The computations of singular values of various matrices (e.g. 
loop transfer, return difference, inverse return difference matrices) were 
carried out using a recently developed multlvarable frequency domain 
analysis software package (FREQ) , and the LQG designs were carried out 
using ORACLS [13]. The nominal plant includes three rotational rigid-body 
modes and the first ten elastic modes. We assume three torque actuators; 
hence, the order of B matrix Is 26x3. Assuming three attitude sensors 
(one for each axis) at the same location as the actuators, C is a 3 x 26 
matrix. The plant, input, and output matrices were obtained from a finite 
element analysis of the antenna. 
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Before starting the controller design, the maximum and minimum 
singular values (o and a) of the full, nominal, open-loop plant transfer 
matrix were obtained and are shown in figure 4. The g plot clearly shows 
the peaks at the elastic mode frequencies (i.e. the poles), the most 

prominent being the first mode near .75 rad/Bec. The dips in g_ indicate 
the presence of transmission zeros for the multivariable plant at those 
frequencies. The controller synthesis studies are performed using the 
design model consisting of: 

a) rigid-body model (n ■ 6, n q * 0) 

b) rigid-body and the first flexible mode (n ** 8, n q M 1) 

c) rigid-body and the first three flexible modes (n - 12, n q “3) 

The measurements available are the three attitude angles at the feed 

location. One 3-axis torque actuator is used at the same location. The 
compensator is designed based on these sensors and actuators. 

5.1 Rigid Body Model : 

In this section the controller design is carried out based only on 
the rigid body design model. The largest and the smallest singular values 
of the rigid-body transfer matrix (n » 6) are of the form l/s^. The 

corresponding additive uncertainty AG, which consists of the (20th order) 
flexible dynamics, is plotted in figure 5. Figure 5 indicates the 

presence of poles at the undamped flexible mode frequencies of 0.75 
rad/sec, 1.35 rad/sec, etc. Also, the pole of the first mode frequency of 
0.75 rad/sec produces the highest peak since it is most lightly damped. 
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(The importance of this fact will be Been later when the stability 
condition is violated at this point). 

For this sixth order design model, a compensator design was carried 
out using the Kalman filter design methodology to achieve satisfactory 
performance (l.e., large gain and bandwidth) at low frequencies, and 

robustness at high frequencies. This design was carried out using the 
Kalman filter Riccati equation (16). The Kalraan-Bucy filter (KBF) 

a 

transfer matrix G KF (s) is given in equation (18). Appropriate 

loop-shaping can be accomplished by proper choice of the weights p and L 
in equation (16). Since the controller design model is of the form 

1/s^, one can analytically evaluate the singular values of I+Gvp 

using equations (16) and (17). Assuming p ■ 1 and L * C^p^)^ the 

left hand side of (15) can be solved. For “ 0 and L£ * ^ 1, it 

can be t'hown that equation (15) ia satisfied by: 

k 2 < 10 ~ 7 

This implies that the Kalman filter gain computed using (17) will be very 
low. Figure 5 shows plot® for condition (15) with two L matrices, with 
* 0 and V -2 " 10”6 and 10“^. The right hand side of (15) is 

also plotted in figure 5. It is evident that condition (15) is satisfied 
for k2“10”7. is decreased further, the curve shifts upward 


thus increasing the margin. 
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The next step consists of LQ regulator design. Having obtained an 

* 

acceptable compensator through Kalman- Bucy filter equations, the LQ 
regulator is realized via the loop transfer recovery method [8). Figure 6 
presents the singular value plots of the complete loop transfer matrix 
GpG c (s) (which consists of the plant, the KBF and the LQR) for 
different weighting parameter q (Eq. 19). The q selected was q-10^ and 
lQ^. It is easy to check condition (11) in this case. As q is 

increased, the plots approach those of the compensator obtained from the 
Kalman filter design approach. The LQ design for q»10^ was considered 
to be satisfactory. 

The standard LQG/LTR procedure requires the definition of the 
"desired" loop transfer characteristics (see step two in section three,) 
That is, o(G KF ) must satisfy the low-frequency performance 

specifications, and a (G^p) must satisfy the high-frequency robustness 
specifications. Thus, in the presence of additive uncertainty AG, the 
procedure states that the robustness condition 


1 d + c K r> Z ( V 

3 <g kf> 


> a (AG) 


should be satisfied. However, in the case described above, it was found 
that the above condition makes the "desired" design (G^p) extremely 
conservative. From f,.«.;ure 6, it is seen that the closed loop bandwith is 
quite low and nowhere near the desired value of .1 rad/sec. Therefore, 


recovering this conservative loop gain yields a compensator with poor 
performance. This fact led to a modification of the LQG/LTR procedure. 
In particular, the above robustness test on Gj^p is omitted in the 
modified procedure. Instead, the recovev j is carried out first, and then 
the (leas conservative) stability test (11) is applied directly for the 
compensator G c » The Kalman filter transfer matrix G^p is based only 
on the desired performance and not on satisfying the stability test of 
equation (15). 

With the revised test on G c > the following choices on L and p 
matrices were made. 


L 



Using the recovery procedure, the compensator is obtained for this 
case wth q « 10^. The resulting stability test (Eq= 11) is shown in 
figure 7. It is seen that the stability margin is lowest at the first 
mode frequency (0.75 rad/sec). Any increase in the gain (obtained by 
q > 10^) resulted in violation of stability condition at that point. 
The overall loop bandwidth 'is obtained from the singular values of the 
loop transfer function G p G c shown in figure 8. It is seen that the 
bandwidth (i.e., the frequency at which cr(G p G c ) = 1) is far short of 
the required 0.1 rad/sec. In order to increase the bandwidth, the gain 
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has to be increased by increasing q. However, this results in the 
violation of the stability condition (11). Thus it is evident that, with 
a rigid-body design model, it is not possible to meet the performance 
specifications. 


5.2 1 Flexible Mode Design Model : 


To overcome the above problems, the next alternative that was 
considered was whether the inclusion of the first flexible mode (0.75 
rad/sec) in the design model would improve the performance. The inclusion 
of the first flexible mode, which is predominantly a torsion mode, results 
in a design model of order 8. The singular values of AG shown in figure 9 
are an order of magnitude lower* than those in figure 5 (wherein AG 
consisted of all the flexible modes). The first peak of o(AG) occurs at 
1.35 rad/sec, which is the frequency of the second mode. This is the 
critical point in the stability test (Eq. 11). After a number of trials, 
the following choice of L and p was made to obtain the desired performance 

G KP* 


lO" 1 !, 


K) -1 x 2 J 0 


The recovery is obtained for q ■ 10 5 and the stability test is shown in 
figure 9. Fig. 9 indicates the critical point to be at about 0.28 
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rad/sec. There ib a good margin at the peakB of AG due to upward sloping 
of the upper curve. The resulting loop transfer function (G p G c ) plots 
are shown in figure 10. The plots indicate that the required 0.1 rad/sec 
/, bandwidth is not obtained (although it is much higher thav. the rigid-model 

case). Any increase in the gain (for q > 10-*) was found to result in 
the violation of the stability condition (11). Figure 10 indicates the 
presence of an open-loop invariant zero near 0.082 rad/sec, which was also 
confirmed by independent computations. This zero is almost on the 
imaginary axis (i.e., the transfer matrix is close to being nonminiraura 
phase). Therefore, (as would be expected) the recovery procedure is not 
very effective for making GpG c approximate G^. 

5.3 3 Flexible Mode Model : 

» 

In order to improve the performance further, the next step was to 

include the first three flexible modes in the design model. It is logical 
to do this because they represent the first modes about each axis, i.e., 
the first torsion mode, and the first bending modes in the XZ and YZ 

planes. Thus, the order of the design model was 12. The singular value 
plots for G p and AG are shown in figures 11 and 12, respectively. It is 
seen from figure 11 that G p has zeros near 0.082 and 0.22 rad/sec, and 
poles near 0.75, 1.35, and 1.7 rad/sec. It is seen from the AG plot 

(Fig. 12) that a is considerably lower than that in figures 7 and 9. 

After numerous trials, the following choice of the L matrix and the scalar 
y was made for a suitable 
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The recovery was accomplished with q - 10*®. The stability test is 

shown in figure 12. It can be seen that condition (11) is satisfied with 
a wide margin. Also, at the peak for AG (at 8 rad/sec) the upper curve 
slopes upward, indicating good tolerance of high-frequency uncertainty. 
The limit for increasing the gain (indicated by the lowest point in the 
upper curve in figure (12) occurs at about 0.3 rad/sec. The resulting 
compensator G c is shown in figure 13. The gain of G c is much higher 

than that obtained in the previous cases. Generally, the LQG/LTR 

technique attempts to choose G c in such a way that the product GpG c 
is replaced by G^p (i.e. G c is attempting to invert G p in the 
frequency range of interest). The 3~mode design plant shown in figure 11 
has elastic mode eigenvalues at -.0075 ± j.75, ■’-.0135 ± jl.35, and -.0170 
± jl.70. Figure 13 shows that G c has zeros with frequencies near these 
locations. The design plant also has transmission zeros at -.9 x 10"^ ± 
j .082 , -.37 x 10~ 3 ± j .22 , and -.29 x 10" 3 ± j.22. Ideally, G c 
should have poles with frequencies near .082 and .22. However, the design 
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plant zeros are too near the jtii-axis and tend to numerically behave as 

» 

nonminimum phase. Some attenuation Is obtained by the compensator pole 

near .4 rad/sec. The plots for the loop transfer matrix G p G c are 
given in figure 14. It is seen that a bandwidth of 0.1 rad/sec. is 

obtained except for the presence of the invariant zero near 0.082 

rad/sec. which causes some deterioration of performance. At frequencies 
past .4 rad/sec. , G p G c behaves like G^p and eventually rolls off at 
60db/decade. Also, a_ and a are closely spaced, indicating good system 
behavior. Thus it is seen that the inclusion of the first three modes in 
the design model yields a robust compensator which also meets the 
bandwidth specifications. 

The final step is to check -the stability of the complete nominal 
system when the compensator G c (s) designed above is used. The overall 

closed-loop system is: 
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where the subscript F is used to denote the full-order nominal plant, and 
x denotes the state estimate for the design model. The eigenvalues of the 
overall closed-loop system using the 3-mode controller are given in Table 
II. It can be seen from Table II that the overall closed-loop system is 
stable. 
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* 

6. CONCLUDING REMARKS 

The LQG/LTR multivariable frequency domain technique was employed in 
the design of an attitude control system for a large flexible space 

antenna. The LQG/LTR method was noted to be especially attractive for 
overcoming spillover effects common to large space structures control 
problem modelled from finite element data. The design objective of 
avoiding excitation of higher order modes while satisfying performance 
criteria was met by including these modes in the robustneos uncertainity 
barrier. 

Design was based on a reduced order model chosen as the rigid body 
dynamics plus the fewest number of low frequency vibrational modes 

necessary to meet a desired closed loop bandwith. Inclusion of the first 
three vibrational modes (corresponding to the three axes) was found to be 
necssary to meet a 0.1 rad/see bandwith. For wider bandwidths, design 
models width greater than three modes may be needed. A satisfactory 

control design was obtained using only a colocated single pair of 3-axis 
attitude sensor and torque actuator for the Hoop/Column antenna problem. 

Performance degradation was observed due to the presence of invariant 
zeros within the design bandwith. These zeros were unavoidable given the 
prescribed sensor/actuator locations and emphasized the fact that 

consideration should be given to control aspects when building large space 
structures. 

A modification of the standard LQG/LTR procedure was introduced in 
which the robustness test was performed with the full LQG compensator 


58 


24 


instead of the intermediate Kalman filter design. This approach was found 
to produce higher gain compensators and helped overcome the basic 
conservativeness shortcoming of the LQG/LTR approach. 
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Table I. Antenna Parameter 8 

Rigid”body parameters 
MasB-4544.3 Kg. 

2 

Inertia about axes through center of mass (Kg-tn ) 

I - 5.724 x 10 6 I - 5.747 x 10 6 
xx yy 

I - 4,383 x 10 6 
zz 

I “I "I "0 
yz xz xy 

Structural Mode Frequencies (rad/sec) 

0.75, 1.35, 1.7, 3.18, 4.53, 5.59, 5.78, 6.84, 7.4, 8.78 
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Table II. 
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Eigenvalues of the Full Closed-Loop system 


Real part 

Imaginary Part 

-8.535 

<1<T 3 ) 

8.054 

do” 2 ) 

-8.535 

do" 3 ) 

-8.054 

do" 2 ) 

-7.557 

do" 2 ) 

1.250 

do" 1 ) 

-7.557 

do" 2 ) 

-1.250 

do” 1 ) 

-7.604 

CIO" 2 ) 

1.248 

do” 1 ) 

-7.604 

Ucf 2 ) 

-1.248 

do" 1 ) 

-2.237 

do" 1 ) 

2.236 

do" 1 ) 

-2.237 

do"* 1 ) 

-2.236 

do" 1 ) 

-2.330 

do** 1 ) 

2.154 

cio" 1 ) 

-2.330 

CIO" 1 ) 

-2.154 

CIO" 1 ) 

-2.379 

CIO" 1 ) 

2.113 

do" 1 ) 

379 

do” 1 ) 

-2.113 

do” 1 ) 

-7.466 

do” 3 ) 

7.466 

do” 1 ) 

-7.466 

do” 3 ) 

-7.466 

do” 1 ) 

-1.346 

Cio” 2 ) 

1.346 


-1.346 

do” 2 ) 

-1.346 


-3.076 

cio” 1 )’ 

1.373 


-3.076 

do" 1 ) 

-1.373 


-1.016 


1.267 


-1.016 


-1.267 


-1.702 

do” 2 ) 

1.702 


-1.702 

cio" 2 ) 

-1.702 


-4.028 

CIO” 1 ) 

1.737 


-4.028 

cio” 1 ) 

-1.737 


-3.181 

Cio” 2 ) 

3.181 


-3.181 

CIO" 2 ) 

-3.181 


-4.422 

Cio” 2 ) 

4.529 


-4.422 

do" 2 ) 

-4.529 


-5.579 

do” 2 ) 

5.590 


-5.579 

cio” 2 ) 

-5.590 


-5.731 

do” 2 ) 

5.776 


-5.731 

do” 2 ) 

-5.776 


-6.685 

do” 2 ) 

6.841 


-6.685 

do" 2 ) 

-6.841 


-6.390 

do” 2 ) 

7.401 


-6.390 

do" 2 ) 

-7.401 


-8 . 326 

do" 2 ) 

8.782 


-8.326 

do" 2 ) 

-8.782 



63 



DEFINITION OF UNCERTAINTIES 


tV 



(b) Additive uncertainity 













Figure 4.- Singular values of the full nominal plant 



STABILITY TEST USING KALMAN FILTER 



Figure 5.- Stability robustness test (equn.15) for KBF-based loop shaping 














Figure 12 . Stability Robustness test (eqn. 11 ) for 3 mode model 






